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Common-envelope events (CEEs), during which two stars temporarily orbit 
within a shared envelope, are believed to be vital for the formation of a wide 
range of close binaries. For decades, the only evidence that CEEs actually 
occur has been indirect, based on the existence of systems that could not be 
otherwise explained. Here we propose a direct observational signature of CEE 
arising from a physical model where emission from matter ejected in a CEE is 
controlled by a recombination front as the matter cools. The natural range of 
timescales and energies from this model, as well the expected colors, light- 
curve shapes, ejection velocities and event rate, match those of a recently- 
recognized class of red transient outbursts. 

Many binary star systems, including X-ray binaries, cataclysmic variables, close double- 
neutron stars, and the potential progenitors of Type la supernovae and short-duration 7-ray 
bursts, are thought to be formed by CEEs. Because most stellar-mass binary merger sources 
for gravitational waves have experienced a CEE in their past, improved knowledge of CEEs 
should decrease the large uncertainty in theoretically-predicted merger rates. However, the 
short timescale expected for CEEs suggested that we would never directly observe them, 
allowing us only to draw inferences from the systems produced. 

A CEE begins when a binary orbit becomes unstable and decays. This might, for example, 
be driven purely by tidal forces (i.e. the Darwin instability), although CEEs are more commonly 
imagined as following a period of rapid mass transfer from one star to the other (1). In some 
cases the rate of transfer is so high that the receiving star is unable to accrete all the matter 
without forming a shared common envelope (CE) around the binary. This CE causes drag 
on one or both stars and hence orbital decay, with orbital energy and angular momentum 
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being transferred to the CE. This may end with a stellar merger or - if the CE is ejected - 
the binary can survive, typically with a much reduced orbital separation, critical to explaining 
many observed compact binaries. 

When a CEE results in formation of a close binary, it is expected that a substantial proportion 
of the mass is ejected - typically almost the entire envelope of one of the stars. Some mass 
can also be ejected in the case of a merger. This partial ejection has two causes (2). Firstly, 
the orbital energy deposited into the CE early in the merger may exceed the binding energy 
of the outer layers. Secondly, the upper layers of the envelope absorb a substantial amount of 
the orbital angular momentum of the companion, and angular momentum transport may be too 
slow to share this across the envelope as a whole. 

Here we consider the behavior of this ejected matter to try to predict the appearance of 
CEEs. A situation involving similar physics -type IIP supernovae-has been previously studied 
(e.g, (3-5)). In that model, as the ejected stellar plasma expands and cools, recombination 
changes its opacity, leading to the propagation of a photosphere-defining "cooling wave," which 
moves inwards with respect to the mass variable. 

For smooth and spherically-symmetric ejecta distributions, the model light-curve will have a 
plateau shape: the area of the photosphere is defined by recombination, and so the emitting 
surface does not grow with the speed at which the ejected matter itself moves. During this 
phase, whilst material ejected by the CEE will expand with velocity of the order of magnitude 
of the initial escape velocity, the photospheric radius should appear almost constant. The 
luminosity Lp of the emission during the plateau (3-5), re-scaled to the likely energy range of 
CEE, is 
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where /?init is the initial radius, is the kinetic energy that the unbound mass munb has at late 
times after escaping the potential well, k is the opacity of the ionized ejecta, and r^c is the 
recombination temperature. The duration of the plateau tp with the same assumptions is 
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This model does not depend on the origin of the energy released during the outburst. For 
type IIP supernovae, recombination controls the release of the internal energy generated by 
strong supernova shocks. For CEEs, however, there is no such supernova-provided energy 
input. Instead the energy released by recombination itself may dominate the energy budget 
of many outbursts (6). The unbound mass munb could potentially radiate - simply due to 
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recombination - as mucin energy as 

Srecom ^2.6x lO'' ei'gS (X + 1 .5F/He) ^ • (3) 

Here X is tine mass fraction of liydrogen and Y is tlie mass fraction of lielium. Hydrogen would 
initially be ionized in almost all of the likely ejected material from most stars; however helium 
may be fully ionized only in some fraction of it, denoted /He- The role of recombination in a CEE 
has hitherto been a debated issue in the overall energy balance, the controversy arising from 
whether it can be effectively converted into mechanical energy to help eject the CE (7-9). This 
energy budget for the outburst may be increased by the thermal energy of the ejecta. Much of 
the pre-CEE thermal energy of the ejecta may be expended on adiabatic cooling (6). However, 
the shock-heating caused by the CEE could well be substantial in some cases. 

We now estimate the extent of the parameter space of CEE outbursts, using the model 
described above to predict the diversity of real events. We assume that scales with the 
gravitational potential at the surface of the primary star (2), and use the dimensionless factor 
^ to write E^ = C,{Gm\fi-a)/Rmiu where /„, = munb/m is the fraction of the total primary mass 
mi that becomes unbound. From equations Q] and [2l this leads to Lp oc {f^mlRrJ^l^^)^^^ and 
?P (/^mi7??„ij(^"')'^*'. Two families of events seem likely, one for mergers (i.e. /m < 1), and 
one for CE ejection (i.e. /m < 1) (Fig. 1). 

In addition to the predicted ranges of outburst energy and duration, this model for CEE 
outbursts has several noteworthy features. The physics that causes the plateau-shaped light- 
curve should lead to a difference in the phometrically-inferred expansion velocity and the actual 
material velocity (which could be inferred from spectra). The effective photospheric temperature 
should be ~ 5000K for thick ejecta (4), and so the outburst color will naturally be red. In 
addition, once the ejected envelope has fully recombined, the material may suddenly become 
transparent — unless enough of the ejecta has cooled down sufficiently to produce dust. These 
characteristics are reminiscent of curious transients with predominantly red spectra recently 
detected in the local universe (e.g., {10-17})). This empirical class has been dubbed Luminous 
Red Novae (LRNe), a subset of the even more ambiguously defined class of intermediate 
luminosity red transients (ILRTs) (2). ILRTs cover a wide range of outburst energies - from 10"^^ 
to a few 10"*^ ergs (brighter than the brightest novae, but still fainter than Type la supernovae). 
They are characterized by spectroscopically inferred expansion velocities of 200-1000 km/s - 
much lower than would be expected for novae or supernovae, and also strikingly different from 
the photometric expansion velocities (18). In addition, some could be seen as red giants within 
a dozen years after the outburst ( 16, 19). 

It was not known what ILRTs are or whether they have a common cause; several ideas 
have been suggested (2). A model which considered the possibility that LRNe are caused 
by stellar mergers - a subset of CEEs - has been independently considered several times 
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Figure 1: Left panel: Model diversity in the Lp-tp parameter space is indicated by lines 
representing constant primary mass and radius, fm is fractional mass loss and C is the kinetic 
energy at infinity, parametrized as a fraction of the binding energy at the surface of the primary 
star. Stellar mergers are in a regime of little mass ejection, while /m = 0.9 approximates full 
envelope ejection. Right panel: Estimated ranges of the plateau luminosity Lp and duration tp 
for primary stars with ZAMS masses from 1 to 150 M©. niunb is the ejecta mass. It is assumed 
that mergers can happen anytime during the primary's evolution, whereas full envelope ejection 
can occur only for post-MS primary stars. We used fitting formulae for stellar evolution (24), at 
Zq. In both panels, values for Lp and tp are marked for the outbursts from V1309 Sco, M85 OT, 
M31-RV, and V838 Mon. 
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for different LRN outbursts, though further examinations of outburst features always showed 
various drawbacks. However, those problematic features do match expectations from our CEE- 
driven outburst model (2). 

A particular feature of the LRN outbursts - as opposed to all ILRTs - is the presence of 
a plateau in their luminosity curves. We compare well-known LRNe (2) to the expected CEE 
diversity in Fig. 1. The agreement is striking, especially given the simplicity of the model and 
the potential complexities it neglects - e.g., how CEE ejecta deviate from spherical symmetry, 
or how much C, for mergers might be different to C, for full envelope ejection (2). 

M85 OT2006-1 is a LRN with well-known peak luminosity and plateau duration. If the 
luminosity from M85 OT2006-1 was largely from recombination, ~ 1.5Mq of plasma would have 
recombined to provide the observed total energy. This fits with constraints on the progenitor 
mass (< 2Mq) from the stellar population age {20). Thus M85 OT2006-1 plausibly ejected the 
whole envelope of a low-mass giant. This outburst showed a plateau, with luminosity ^ 1- 
2 X IO^Lq {21) for sa60-80 days, and had expansion velocities of 350 km/s {14). Our inferred 
ejecta mass and that observed expansion velocity indicate a kinetic energy of ~ 1.8 x 10"^^ erg. 
Then Ri^n = 45/?©, self-consistent with our model, gives Lp ~ lO^Lg and tp ~ 70 days. 

Another recent outburst, V1 309 Sco, is similar to, but fainter than, most LRNe, as it radiated 
away only ~ 3 x 10"^"^ erg during a ~ 25-day plateau-shaped maximum in the light-curve {19). 
The progenitor was a contact binary with a relatively rapidly-decaying orbital period of ~ 1.4 
day. After the outburst, the system appeared to be a single star; therefore this appears to have 
been a CEE, leading to a merger (79). However, several features of the V1309 Sco outburst - 
in particular the plateau in the light-curve and sudden transparency - were difficult to reconcile 
with prior theoretical expectations for the appearance of a CEE (2). 

Because the V1309 Sco progenitor was observed in detail, this system is ideal for testing our 
model. Beginning with the properties of the pre-merger contact binary {19,22), we calculated 
the amount of material that became unbound during the V1309 merger using two methods - 
simple energy balance using a 1 D stellar code and a set of 3D hydrodynamical simulations (2). 
Both methods predict that a small mass, ~0.03 to O.OSMq, will become unbound. Complete 
recombination of this ejected mass would provide enough energy (> 7 x 10"^"* ergs) to explain the 
total energy output of V1309 Sco. The output of the 3D simulations, combined with Equations 
[Hand 121 predicts plateau durations from 16 to 25 days, and plateau luminosities of 1.8-2.2 x 
IO'^Lq. These values match the observed luminosity (Lp ~ 2 x IO^Lq) and plateau duration 
(about 25 days). 

Considering the menagerie of theoretically-expected outbursts from CEEs, we note that 
events in the top-right of Fig. 1 should be relatively rare [compare with rj Car; see (2)], and those 
in the bottom-left (stellar mergers) comparatively hard to detect in a magnitude-limited survey. 
Assuming that the peak luminosity of the outburst is about an order of magnitude higher than 
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Lp, we find that the whole range of Lp and tp for stellar masses 1 - 15OM0 coincides well with the 
observed domain for luminosities and durations of LRNe suggested in (23). We can estimate 
the rate of CEE-originated outbursts that appear as red transients, by considering what fraction 
of stars in the galaxy undergo a CEE. We estimate 0.024 such events per year per Milky Way- 
like galaxy (2), of which about half should be more luminous outbursts (results of a CE ejection), 
and half are lower-luminosity events (powered by stellar mergers). This is consistent with the 
empirical lower limit for more luminous ILRTs of 0.019 yr^ for the Galaxy (20), because we do 
not expect that all luminous ILRTs must be powered by a CEE [though some non-LRN ILRTs- 
like NGC 300-OT or SN2008S - might potentially also be triggered by CEEs (2)]. 

The question of whether recombination energy can help to unbind a stellar envelope during 
a CEE is important for understanding the formation and survival of many binary systems {8, 9). 
Our model suggests that a large fraction of the energy from recombination is commonly radiated 
away following a CEE. Such luminosity provides a beacon, which helps to illuminate and identify 
a CE ejection or merger at large distances. The recombination wave also controls the shape of 
the plateau-shaped light-curve of LRNe. We therefore suggest that detecting and characterizing 
the population of ILRTs will help us understand CEEs. 

Supplementary text 

Here we provide additional information about our hydrodynamical simulations, as well as 
give more details behind our formulae and arguments from the main text. 

S1. Intermediate Luminosity Red Transients 

In recent years, observers have identified a new class of transients, with peak luminosities 
somewhere between that of the brightest novae and Type la SNe and with a total energy output 
anywhere from 10"^^ to a few 10"^^ ergs, e.g., (10-17). Because their spectra are predominantly 
red, completely unlike novae, they have been called intermediate luminosity red transients 
(ILRTs, (25)). However alternate names are also in common use, either for all of the transients 
or for particular kinds of events, e.g. Luminous Red Novae (LRNe), intermediate luminosity 
optical transients, mid infra-red transients, V838 Mon-like events, and supernova impostors. 
The relative novelty of these classes, and the uncertainty over their physical origins, makes it 
unsurprising that the use of such terms is not necessarily always consistent between different 
papers. 

It is not known what red transients are or whether they have a common cause: it has 
been argued that ILRTs could be due to accretion-powered jets (26), tidal disruption of planets 
(27), nuclear outbursts (28), electron-capture SNe (29) with shocks propagating through dusty 
surroundings (30), or violent stellar mergers (31). It has also been argued that - despite having 
similar outburst properties in terms of color and luminosity - ILRTs may have different origins 
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(25). 

In this work, we consider the particular group of red transients which are most frequently 
labeled as LRNe. The widely-recognized members of this class of red transients are M31-RV, 
V838 Mon, and M85 OT2006-1. We note that the recognition of LRNe as a special class of 
stellar explosion started with the last object in that list, although it was not the first of them to be 
detected. In this paper we also include V1 309 Sco as a possible member of this LRN class, at 
the low-energy end ( 19). Another suggested low-energy end LRN is V4332 Sgr - this transient 
radiated away a total energy of only 4.5 x 10"^^ ergs (32), an order of magnitude less even 
than VI 309 Sco. The LRN class is likely separated from another class of red transients where 
progenitors were observationally identified to be dusty modestly-massive stars (with M lOM©) 
- as for example in the cases of SN 2008S and NGC 300 OT ( 18,33,34); those ILRTs are better 
known as supernova impostors. Nonetheless, because the recognition of red transients as a 
new type of astrophysical object is very recent, their classification is not yet by any means 
universally established or finally accepted throughout the astrophysics community. 

We note that some investigators have suggested that LRNe could be present only in old 
populations, with V838 Mon standing as an exception (see also the discussion against this 
point of view in (35)). It is important to clarify that our proposed link between CEEs and LRNe 
implies that there should not be any such restriction on population age for LRNe: our CEE 
outburst model is related only to the reason of the ejection - a CEE - and to the self-similar 
physics and recombination outburst following the CEE. Indeed, CEEs may occur in binary stars 
with companions of any masses. In terms of rates, CEE from massive stars are expected to be 
much more rare, as their relative fraction in the total number of stars is small; however they are 
not forbidden. In addition, given the current understanding of CEE physics, it is impossible to 
say how the CE ejection would differ in low-mass and high-mass stars. 

S2. Previous problems with a merger model: what does the standard CEE model 
predict? 

One of the most common models suggested to explain different LRNe - V838 Mon, M85 
OT 2006-1 and V1309 Sco - is a merger, either of two stars, or of a star with a planet (e.g., 
{13, 19,31,39)). It is also known as the merger-burst model; however, because a merger is a 
particular type of CEE, we will hereafter refer to it as to a CEE model. The case of VI 309 Sco 
is rather unequivocally a CEE, as observations revealed a binary prior the outburst and a single 
star after (19). However, there were problems with justifying the physics behind this model. 
For example, a stellar merger scenario in the case of M85 OT2006-1 was rejected because an 
analysis of the energy available from a merger suggested that a violent merger would be unable, 
by a factor of a few, to explain the energy that was radiated away from M85 OT2006-1 {20). 
In (37), striking similarities between the M31-RV, V838 Mon and V4332 Sgr outbursts were 
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noticed, leading to a discussion wliicli concluded that the merger-powered outburst models 
proposed to-that-date showed too much dependency on metallicty, mass, and ages and hence 
could not explain the observed homology. 

We list next a number of striking observational features of LRNe which were not previously 
theoretically anticipated by a CEE model (one leading to either a merger or a the formation of 
a close binary) and which were not emphasized as discriminating features in previous studies: 

• Large increases in radius and iuminosity. In previously published CE simulations, 
the increase in the stellar radius of the bound mass during fast spiral-in is usually less 
than a factor of 10 (see e.g., (40,41)). These simulations were performed for binaries 
that are expected to survive, rather than merge, i.e. these binaries have more energy 
injected into the stellar envelope - relative to the binding energy - than for binaries that 
merge. However the well-studied example V 1309 Sco reached a maximum luminosity 
~ 3 X IO'^Lq, about ^ 5000 times larger than the initial luminosity of the progenitor (~ 3.0- 
8.6L0), whilst its effective temperature dropped to ~ 4000K during the peak luminosity 
and plateau {19). This indicates an increase in its apparent (or effective) radius by about 
90 times, reaching ~ 300Rq (19), at least an order of magnitude more than simulations 
suggested. Similarly, interferometric observation showed that V838 Mon expanded by 
^ 300 times during the outburst, to 1570±400/?q, while its progenitor was likely 5 - IOM0 
with a radius of only about 5Rq {42). 

• An unexpectediy iong duration. Expansion of the envelope during CEE is expected to 
occur during the dynamical plunge-in phase (see, e.g., (1)), i.e. the stage during which 
the inspiralling star loses most of its angular momentum. It is theoretically expected, 
and confirmed by numerical simulations (see e.g. {40,41)), that this phase proceeds on 
a timescale of about several initial dynamical timescales of the stellar envelope of the 
primary star. This is a few days for the case of VI 309 Sco; however, the outburst lasts 
sal 00 days ( (19), also see Fig. SI) which is comparable to dynamical timescale of this 
object at the maximum of expansion, Tdyn.max ~ 80 days. In the case of V838 Mon, the 
outburst last for about 80 days (43) while the initial dynamical timescale would be only a 
couple of hours. In case of M85 OT2006-1 , in which the outburst last about 60 days ( 14), 
there are no strong constraints on progenitor. 

• Piateau phiase. Several ILRTs outbursts, after the initial rise in brightness, featured a 
plateau in their light curves in the red band at a luminosity somewhat lower than at the 
maximum. Specifically, in case of VI 309 Sco the plateau duration is ~ 25 days (depending 
on the definition of the start and the end of the plateau phase, from 1 6 to 31 days ( ( 19)), 
also see Fig. SI), -60-80 days in case of M85 OT2006-1 (14), ~ 15 days in case of 
V4332 Sgr ( (32), although the start of the outburst is not well known) and ~ 30-40 days in 
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case of M31-RV (37). V838 Mon has a very complex visual band light-curve with several 
peaks separated by about 30 days. Nonethless, it does have an easily distinguished 
plateau in the red part of the spectrum as well as in bolometric luminosity; the overall 
duration of the outburst is about 80 days, and that plateau starts ~ 30 days after the peak 
luminosity and lasts for another ~ 35 days {36, 44). Most of the energy in LRNe is radiated 
away during the plateau phase. The plateau shape had not been predicted and had not 
yet been explained. 

• An extremely rapid decline. The post-common-envelope decline in luminosity of the 
merger product is expected to occur on a thermal timescale, which - for these over- 
luminous and over-inflated objects - is typically about a year. However, the abrupt decline 
in luminosity (by a factor of 100) in VI 309 Sco, M31-RV, V4332 Sgr and V838 Mon 
happened in a dozen days ( (19,32,36,37,44), also see Fig. SI). For V 1309 Sco, the 
timescale for exponential decrease was determined to be only a few days ( 19). Strikingly, 
the timescale for this sharp luminosity decrease is considerably shorter than the dynamical 
timescale Tdyn,max at the maximum of expansion for V 1 309 Sco (80 days) or for V838 Mon 
(about 500 days). 

• Inconsistent velocities. While expansion velocities from the spectra were several hundred 
km/s (e.g., in VI 309 Sco {19), in M85 OT2006-1 {14), and in V838Mon (45)), velocities 
inferred from the apparent radius expansion imply speeds from only about 20 km/s (V 
1309 Sco) to about 100 km/s (V838 Mon). 

In section S3 below we explain how the model presented in this paper naturally accounts 
for the above features. 

S3. Wavefront of cooling and recombination in a CEE 

As discussed above, the light curves of LRNs have five similar striking features, where two 
are the most exceptional in providing important clues to the physics of the outburst. First, the 
outbursts have a roughly constant plateau luminosity Lp for a time ?p that is typically dozens of 
days (for a note on the special case of V838 Mon see, see §S5). Second, the inferred radius 
of the photosphere increases relatively slowly during most of the plateau phase, at the end 
of which the apparent radius decreases very quickly. This is reminiscent of the behaviour of 
type IIP supernovae, in which the photosphere does not stay at a fixed Lagrangian coordinate 
but moves inwards in mass as the ejecta expands and cools; roughly self-similar, homologous 
expansion means that the radius of the photosphere and luminosity of the emission remain 
approximately constant (4, 46^9). 

Such a photosphere-defining 'cooling wave,' propagating inward in the frame of the expanding 
shock, was predicted for terrestrial explosions in air by {50,51). This type of behaviour does not 
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need to be associated with recombination; it is only dependent on a non-linear temperature 
dependence of the opacity of the ejecta (see also (52), which argues that the necessary 
condition is that the effective opacity increases more rapidly than a cubic power of temperature). 
The low-opacity material outside the photosphere is able to cool by radiation much more 
effectively than the hot, high-opacity material inside. However, in the case considered here, 
as well as for type IIP SNe, recombination is the cause of the change in opacity that leads to 
this cooling-wave. The location of the effective photosphere is controlled by a sudden large 
reduction in the mean opacity of the ejecta after hydrogen recombination. 

Even though the mean opacity greatly drops outside the region of recombination, the photons 
emitted during recombination will not in general escape directly. Although neutral hydrogen 
has a much lower Rosseland mean opacity when compared to ionised hydrogen, the relevant 
line opacity of neutral hydrogen is higher than for the ionised material. So the recombination 
photons themselves have a short mean free path, and there is no reason to expect strong Ha 
line emission, even if photons at other wavelengths are free to escape. Uniformly applying 
the Rosseland mean opacity in this case would be misleading, and full wavelength-dependent 
radiative transfer would be necessary to simulate precisely the structure of the cooling wave. 
Hence situations where the location of the cooling wave is defined by recombination are more 
complicated than the cases first considered by (50,51), in particular because there is a release 
of energy from recombination that is coincident with the change in opacity. However, we 
consider that the energy input is so close to the photosphere that it justifies neglecting the 
effect of any heating on the main features of the hydrodynamics. 

This physics is well-studied in the context of supernovae, and we can apply appropriate 
results here. The analytic model of (4), based on the idea of there being a wavefront of 
cooling and recombination (WCR), is particularly attractive. Although the method produces 
only approximate light curves, it provides accurate scaling relations that give insight into the 
underlying physics and to the interpretation of observations. Indeed, (5) verifies the major 
results of the Popov model with detailed hydrodynamical simulations that include radiative 
transfer. 

The analytic WCR model has the following intrinsic features: 

a) An outburst has a plateau phase. (For more on the condition which must be satisfied in 
order to obtain a plateau, see below.); 

b) Because the observable surface during the outburst is not the stellar surface, but instead 
the photosphere of the expanding ejected material defined by the location of the WCR, a 
large apparent radius increase is expected; 

c) As this recombination front propagates inwards through the material by mass, its geometrical 
distance to the center can gradually increase and then decrease until the recombination 
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of the material is over (see Eq. (16) in (4)), wliile at tine same time tlie outflow passing 
through the WCR continues to stream outward. This leads to what could be observationally 
identified as an inconsistency between measured material expansion velocities and apparent 
expansion of the radiating surface: the photometric velocity indicating the expansion rate 
of the photosphere can be much less than spectroscopic velocities, which give the speed 
of the gas passing through the photosphere; 

d) The model connects the outburst not with the dynamical timescale on which the CE rapidly 
expands, but with the timescale over which recombination wave proceeds; 

e) After the recombination is completed, the model allows sudden transparency, which could 
be observed as rapid decline of luminosity and appearance of the central object. If 
material above the front is dense and cool enough, it could form dust that would hide 
the central object again (for a discussion of dust formation in the ejecta of a related class 
of transients then see, e.g., (53)). 

The recombination process depends strongly on how not only temperature, but also density, 
evolves. Hydrodynamical simulations show that recombination temperature in this model is ~ 
5000-6000K for thick and dense ejected envelopes, as in case of SN IIP (4,5). Somewhat colder 
recombination temperatures are expected in less dense ejecta, where the number density 
of free electrons (which are required for recombination) is smaller. In the grey atmosphere 
approximation, the effective temperature is a factor of 2^/"^ 1.2 larger than the recombination 
temperature, corresponding to effective temperatures in the '--4000- 6000K range (see also (54) 
who showed that effective temperatures shift to lower values once the energy of the outburst is 
decreasing). Although full 3D hydrodynamical simulations with radiative transfer are necessary 
for more detailed calculations of the effective temperature, it is clear that the WCR model 
predicts that a CEE outburst will be bright in the red part of the spectrum. The above is a 
simplification of the complicated situation occurring in these photospheres (e.g. it neglects the 
role of H ions formed due to electrons donated from metals), but the outcome of a detailed 
analysis should yield similar conclusions. 

For this study, we take the WCR model derived by (4), in the convenient form expressed 
by (49), and re-scale it to the situation we are studyingllJ We scaled opacities to 0.32 cm^g^ 
for X = 0.7. The original derivation by (4) assumed that the kinetic energy input produces self- 
similar expansion of an envelope moving with constant velocity. This implies that the potential 
well of the exploding star has already been overcome, and hence that energy input for his 
equations in the case of CEE is the same quantity as E'^. 

''For equation (1) in the main text, we have corrected a small numerical error in the coefficient of equation (9) 
of (49), for the plateau luminosity, such that it is consistent with the (4) value. 
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An important assumption in Popov's model is tliat recombination, altliougli controlling the 
energy release through its effects on opacity, does not dominate the dynamics of the expansion. 
In the case of CEE, we can estimate from Eq. (3) of the main paper how much energy is 
available from recombination. For typical cases of partial envelope ejection, such as in the 
simulated case of V1309 Sco below, the kinetic energy in the outburst is 3-4 x lO'^^ergmunb/^©, 
exceeding the recombination energy by more than an order of magnitude. For cases of complete 
CE ejection, we can estimate the kinetic energy to be of the same order as the initial 
potential energy of the CEE progenitor-star: = Gmunh{mi -n^unb)(A/?init)"^ where A is an 
envelope structure parameter that is roughly of order unity for small low-mass giants but that 
can be as small as 0.01 for massive and larger giants (55). Then 



150x(X + 1.5F/He 



-Erecomb -^0 V ^0 

So E^ easily exceeds i'recomb by at least an order of magnitude. Hence, although the recombination 
energy is a larger fraction of in the CEE problem than it is in the Type IIP SNe problem, 
the recombination energy is still relatively small enough that we feel justified in assuming that 
it should not play a significant role in driving the dynamics, therefore making Popov's model 
applicable for CEE. 

We briefly clarify here that the analysis for E^ above must self-consistently extend to the 
choice of C, for complete envelope ejection. For that reason, a typical C, for complete envelope 
ejection might be expected to be somewhat higher than a typical C, for ejection resulting from 
a merger event. In those extreme cases where A = 0.01 then an assumption of <^ = 100 would 
become reasonable, though the range 1<C, <IQ seems likely to be more typical. 

In Popov's simplified analytical model for the WCR, the main role of recombination is to 
provide the reason that the gas switches from opaque to transparent, resulting in a plateau- 
shaped light-curve. The Popov model predicts that light curves will be self-similar to each 
other, with their shape characterized by the dimensionless parameter A: larger values of A 
correspond to more pronounced plateaus. We find that, scaling for the case of CEE-outbursts, 
A from (4) can be calculated as 

\3.5RqJ yiO^^ergJ \0.03MajJ V0-32cm2g-iy \^ 4500x7 ' 
which for cases of complete envelope ejection becomes 

A ~ 1700AV2 ^^V'^' ( """^ ^ ( ^ V 

^ ''''^ \ReJ \{mi-m^Me) V0.32cmVV WsOOK^ " 

Consequently, for most of the range of possible CEE-progenitors, A > 1, and hence a 
plateau phase with sudden transparency is predicted. If A < 1, which is not a likely case for 



12 



CE progenitors, the model predicts slow decline with no well defined plateau in the light-curve. 
Cases with intermediate values of A ~ 1 may provide a phase somewhat like a plateau, with the 
luminosity taking a comparatively long time to decline after recombination is completed (see 
also Fig. S2). 

We note that these scaling relations were derived under the approximation that the ejection 
is spherically symmetric, which may well be less accurate here than for supernovae. In addition, 
(4) assumed that the plasma was radiation dominated, which may well be less appropriate for 
matter ejected from a stellar merger than for supernova ejecta. However, it has previously 
been noted that at the location of the cooling wave recombination itself will automatically lead 
to radiation becoming dominant (e.g., (3)), which helps to justify our adoption of the Popov 
expressions. Indeed, (6) show that during the recombination phase the third generalized 
adiabatic index (F^) decreases below the normal value for radiation-pressure dominated matter, 
even for initially gas-pressure dominated matter. This allows continued expansion by an order of 
magnitude or more at near-constant temperature - independent of whether the initial matter is 
gas- or radiation-pressure dominated - and also increases the relative importance of radiation 
pressure, even in initially cold envelopes. 

S4. Simulations of V1309 Sco 

For the properties of the pre-merger contact binary in VI 309 Sco, we take a low-mass 
subgiant with mass Mi ~ 1.52Mq and radius Ri = 3.5Rq with a lower-mass companion of M2 ~ 
O.I6M0 (22), with an orbital period at Roche Lobe overflow of Poa ~ 1-42 days (the last detected 
orbital period before the merger, {19)). 

The binding energy of the envelope of a subgiant representing VI 309 Sco primary is \Ebmd I - 
1.5 X 10"*^ ergs (as is standard, here we include both the gravitational potential and internal 
thermal energy terms but not the recombination term). Stellar model calculations were performed 
using the stars /ev stellar evolution code, originally developed by Eggleton (56-59), with 
recent updates described in {60, 61) and references therein. For this companion mass of 
O.I6M0, even if the companion survives the entire spiral-in and the orbital energy is released 
as close as possible to the boundary of the primary's core (i.e. providing the maximum possible 
energy output), the merger would produce only about A^oa — 5 x lO'*^ erg. This is only a third of 
the envelope's binding energy. So a CEE in this system necessarily results in a merger before 
the envelope is unbound. 

However, for this particular primary, soon after the companion starts to plunge-in, the 
released orbital energy near the surface exceeds the binding energy of that part of envelope 
already outside the location of the companion. From the detailed stellar model, we find that the 
mass that can become unbound during this initial phase of the spiral-in is wiunb ~ O.O4M0. The 
total orbital energy that was deposited in this mass is ~ 6 x 10'*'' ergs; some of this material will 
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be ejected only barely faster than the local escape velocity (vesc ~ 420km/s for the unperturbed 
star) and some will get significantly more specific kinetic energy. Ejected matter that is given 
more energy than is required for overcoming the potential barrier will still have non-zero kinetic 
energy at infinity > 0; this will be smaller than the total energy deposited, but should be of 
the same order. 

The magnitude of the initial velocity v of the ejected upper layers can be qualitatively understood 
by considering a circular binary consisting of two companions with masses mi and m2, orbital 
separation a, and total orbital angular momentum / approximated by Keplerian two-body expressions: 
J = 1x0.0^ = mim2(Ga/M)^/^, where the reduced mass = mim2/M, D.^ = GM/a^, and M = 
mi +1112. To shrink the orbit by da, a fraction of the orbital angular momentum dJ must be 
transferred to a portion dmi of the stellar envelope material. This angular momentum transfer 
occurs in the very top layers of the primary star, and hence the stellar mass which is inside 
the orbit of the companion star does not significantly change. Therefore, treating mi and m2 as 
constant. 



When the companion already orbits inside the envelope, a « n (although the n describing 
particular particles of matter can be somewhat larger than their initial ri in an unperturbed 
star). In the process, the stellar material dmi achieves an angular momentum dJ = vridmi, 
implying a tangential velocity 



For the masses and initial orbital period used in our VI 309 Sco simulations, we find that 
this initial tangential velocity is 



From the mass profile of the larger star's envelope, dmi/dri = 47irjp{ri), where p(ri) is the 
density profile and ri is the distance to the center of the donor. In the upper O.O2M0 of our 
pre-merger 1.52 Mq star {dmi/MQ.)/{dri/RQ.) ~ lO^'^-lO"^ with smaller values closer to the 
surface. We note that when mass ejection begins, the layers will expand and somewhat more 
mass will have even smaller values of {dmi/MQ)/{dri/RQ). A significant fraction of this upper 
envelope therefore is expected to be ejected at velocities higher than Vesc- Considering the 
pre-merger orbital configuration, the expected kinetic energy at infinity is of order ~ 10"^^ ergs. 

The estimates above were made for a symmetric, non-rotating 1 D stellar model. For a real 
-3D - star filling its Roche lobe, most of the companion-envelope interactions occur close to 
the orbital plane. There the star's radius would be larger than for 1 D models, so we expect that 



dJ ~ ^—da = —ixQ.ada . 



la 2 
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more mass could be lost in the rotating case. Yet, since most of this expansion of the stellar 
structure happens in the very outer layers, we do not expect the enhancement to be more than 
a factor of a few. 

From the stellar model, we also find that fue ~ 0.8. Hence, if all mass that could become 
unbound would also recombine, the total energy that can be radiated away during the recombination 
in V1309 Sco case is ~ 10"^^ ergs - so the available recombination energy reservoir can explain 
the observed outburst's energy very well. 

For a better estimate about how much mass could be ejected, we performed several numerical 
simulations using the 3D SPH code starCrash {62, 63). This code was specifically re- 
developed to deal with close binary systems (64). For these numerical studies, we varied 
the initial orbital period around the observed pre-merger value Porb ~ 1-42 day (19). Table S1 
gives the complete list of initial conditions, including rotational synchronization of the giant. The 
companion was modelled as either a low-mass main sequence star (with SPH particles) or as 
a point mass (representing the core of a red giant which lost its envelope in a previous binary 
interaction (22)). Stellar structures were first calculated using the stars /ev code and then 
relaxed in starCrash in a binary configuration close to Roche-lobe overflow. As a result of 
this relaxation, the radius of the giant in the orbital plane was slightly larger than its 1 D radius 
obtained with the stellar code, as expected (65). 

A visualization of simulation ps334 presents the evolution of the column density (in g/cm^) 
as viewed from a direction perpendicular to the orbital plane (see Movie SI ). After the merger, 
~ O.O6M0 of unbound material is left streaming from the merger product with some deviation 
from axisymmetry, including dumpiness. In all simulations the merger ejects a small fraction 
of the donor, and the unbound mass is comparable with our preliminary estimates, ~ O.O4M0 
(depending on the degree of corotation, it varies from ~ 0.03 to O.OSMq, see Table S1). This 
mass is obtained by computing the total energy (kinetic, gravitational and internal) for each 
particle: if this energy is positive, the particle is considered to be unbound. We find that for the 
same donors, having a MS companion will result in less mass ejection. For most simulations, 
the mass is ejected in two episodes (see Fig. SI ), where during the first episode the companion 
is still outside the giant, and during the second it is plunging-in. In simulations mn351 and 
ms376, three episodes of mass ejection were observed. In one simulation, ms372, the mass 
ejection was mostly continuous. The mass ratio between first and second mass ejection varies 
from about 80% of all unbound mass being ejected during the first mass outburst to about 36%. 
The time interval between the starts of the first and last mass outbursts varies from 1 to 9 days, 
and the specific kinetic energy during the first mass outbursts is larger (see, e.g.. Fig S3). 

As reported in Table S1, the kinetic energy of the unbound material at infinity is ^ 10"^^ 
ergs. The initial energy powering the ejection Ey, is larger, and is partly spent on overcoming 
the potential well just near the merged star, so E^ is not identical to E'k- In each individual 
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simulation, asymptotically approaches some value soon after the merger (see Fig. S3). We 
find that the total value of is larger when the unbound mass is larger, but the specific kinetic 
energy is lower. 

We note that the observed V1309 Sco light-curve (see Fig. S1) might well be reproduced 
best assuming that, as in most simulations, there were two or three episodes of mass ejection, 
where in the first episode the unbound mass had higher ejection velocity (Fig S4). It is tempting 
to interpret this as a promising sign for this model, and not just to the extent that adding 
more free parameters will always allow a better fit. Nonetheless, we recognise that a strong 
conclusion with regard to this point requires more sophisticated radiative transfer calculations 
in order to produce light-curves that properly take into account the asymmetry and structure of 
the ejected matter. For the case of ps334, the asymmetry of the photospheric surface can be 
seen in Movie S2. 

S5. Individual objects 
V838 Mon 

V838Mon's light curve is an exception as being more irregular - its V magnitude exhibits 
three phases of brightening (44), with plateaus being very pronounced in bolometric luminosity 
(which is dominated by red spectra) during the last two peaks in V. This complex behavior 
might easily be linked to a dramatic and asymmetric start of a CEE. One way a GEE event 
can start is because of the Darwin instability (which is likely the case for V1 309 Sco, as it is a 
system presumably with a large mass ratio); in this case a CE is formed almost immediately, 
on a dynamical timescale. The other trigger for a CEE event is when one of the stars starts 
overfilling its Roche lobe due to evolutionary expansion (this might be the case for V838 Mon). 
In the latter case, the creation of a CE does not necessarily occur on a dynamical timescale: at 
the very beginning of the contact, loss of surface layers is not expected to occur smoothly and 
continuously. Mass transfer can bring a giant out of contact, either due to a binary expansion 
or due to contraction of the giant in response to sudden mass loss, and the mass transfer 
unavoidably will be continued later. This stage of evolution - the start of the CEE in giants - 
is not well understood and is currently under debate (for a review on the current understanding 
see (1)). The binary may then proceed to a CEE on a timescale which can be estimated to be 
anywhere from a year to hundreds of years ( 1, 66). We suggest hence that the observed three 
peaks of brightening are most likely linked to either one highly asymmetric interaction resulting 
in many clumps or to several mass outbursts, similar to the case VI 309 Sco discussed above, 
but possibly not resulting in a merger. 

It is remarkable as well that it is likely that V838 Mon is a binary star now (67), and hence 
may deliver several more interesting eruptive events before it completes its CEE evolution with 
either merger or the formation of a compact binary. 
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ILRTs from massive dusty stars 

Some red transients have progenitors well identified to be cool red giants surrounded by 
dust. Noticeable examples are SN 2008S and NGC 300 OT (18,33,34). A very good model 
explaining the physics of an explosion propagating through a cool expanded envelope was 
developed by (30). A GEE outburst model as described in our manuscript might seem - at first 
sight - not applicable. To start with, GEE outbursts presume that an envelope is still hot enough 
to be ionized. 

However, the reason for the explosion in that model (30) is not fully established. It has 
been argued that it could be an electron-capture supernova (18) which occurred in an extreme 
asymptotic giant branch (AGB) star. The progenitor is suggested to have been a cold star with 
Ten = 2500/: and Igio(^/^0) = 4.6 and 4.9, where the total energies radiated away in the transient 
stages were 3 x 10^^ andS x 10^^ erg, for SN 2008S and NGG 300-OT respectively (30). 

Whilst a progenitor with those properties is consistent with being a massive AGB star, similar 
characteristics could also be possessed by a red giant during a long-term quasi-stable GEE 
phase known as the self-regulating spiral-in. This situation is expected to occur if the envelope 
was not ejected promptly during the plunge-in phase; in this case the companion orbits inside 
the GE in a very rarefied region and the phase can last hundreds of years {66). The stars 
appear to be puffed-up and somewhat cooler than would be predicted by standard evolutionary 
tracks (68). In addition, the initial mass outburst during the initial interaction and the plunge-in 
(see the discussion on V838 Mon above) will cool down, potentially forming dust around the 
system. In this case, the main orbital energy release occurs deep inside the star, at the bottom 
of the expanded envelope. For stars orbiting AGB cores, the orbital energy release can easily 
be as high as a few times 10"^^ ergs. 

We therefore note that we can not rule out GEE from the list of possible triggers of the 
explosion underlying this subclass of ILRTs, even though the physics of the outburst is different 
from the WGR-determined model, since the outer parts of the extended envelope are cool 
enough to have recombined before the outburst. 

T7 Carinae 

The Great Eruption of rj Garinae, which may also have been triggered by a stellar merger 
(see, e.g., (69) and references therein), was recently measured via light echoes, and was 
susprisingly cool (70), as would have been produced if a recombination wave defined the 
photosphere. The amount of energy radiated was certainly too large to have been powered 
purely by the recombination of a sensible amount of material, so additional energy input would 
have been required. 
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S6. Rates 

A lower limit for the ILRT rate from observations is 0.019 yr^ L^^ {20), where Lmw 
represents the luminosity of our Galaxy. With a star formation rate of « 2Mq per year (though 
rates can be from 1 to 10 M© per year, (77)) and with initial mass function from (72), we find 
that, per year, roughly 0.3 stars are formed with masses large enough (> IMq) to evolve off 
the MS in less than a Hubble time. Roughly half of stars are in binaries, and, for these initial 
masses, 1 6% of those binaries evolve via CEE (with 48% of the binaries surviving the CEE and 
the other 52% experiencing a merger) (73). Hence the theoretically expected rate of CEEs is 
~ 0.024yr"i L^^, in agreement with the observed lower limit. 
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Fig. SI. Observed evolution of magnitudes or luminosity for LRNe discussed in the text. Data 
was taken from {14) for M85 OT2006-1, from (36) for V838 Mon and from (37) for M31-RV. 
For V1309 Sco, we used data from AAVSVO (available at h ftp ://www. aavso . o rg/j , marked by 
asterisk symbols, and from OGLE-II (38), marked as solid dots. 
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Fig S2. Values of A, i.e. the parameter which indicates whether the outburst is expected to 
lead to a plateau, as calculated from (4). We show how these A-values vary across a notional 
pre-CEE mass-radius plane for a selection of parameter choices indicative of different potential 
physical situations. The two left-hand panels are estimates for mergers which eject one percent 
of the mass of the primary star (i.e. /m = 0.01); the lower panel adopts ^ = 0.3 and the upper 
panel ^ = 3. The two right-hand panels represent estimates for envelope ejection (i.e. /m = 
0.9); the lower panel adopts C = 1- and the upper panel ^ = 10 (as discussed in the text, it is 
reasonable to expect that may typically be higher for full envelope ejection than for mergers). 
Only when small amounts of matter are lost from unusually large low-mass stars does A fall 
as low as 1 . Hence we expect a plateau to be produced in the vast majority of CEE-related 
outburst light-curves. 

20 




Fig. S3. The time evolution of tine kinetic energy of tlie ejected material (top panel) and 
unbound mass munb (bottom panel) in simulation ps379. In the lower panel, note the two 
episodes of mass ejection, near times of 6 days and 10 days. 
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Fig. S4. The observed evolution of bolometric luminosity during the VI 309 Sco outburst ( 19) is 
represented by solid circles and a thin line. We compare this to the solutions of equations (17) 
and (21) of (4), using two consecutive mass outbursts. The light-curve resulting from the first 
mass ejection is shown by squares: munb = O.O2M0, = 0.9 x 10"^^ ergs, Ej^ = 0.13 x lO'^^ is the 
initial thermal energy in this ejected layer, Eth = Ek/2+Ef}^ (the material is shocked, the standard 
case for Popov's model when Ej^ < ^k), the fit was done with rrec=5000 K, characteristic for 
this peak (79). The light-curve resulting from the second mass ejection is shown with triangles 
munb = O.O4M0, £k = 0.9 x 10"^^ ergs, = 0-75 x 10"^^ is the initial thermal energy in this ejected 
layer, = shock heating in second ejection. We adopted this for two reasons: a) 

this layer possesses Ej^ larger than it would be suggested by Popov's approximation £'k/2; b) 
the second mass outburst occurs after the companion is below the ejected layer), the fit used 
7;ec=4000 K, characteristic for this peak (19). The interval between the beginnings of the two 
outbursts is 6.7 days. 
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Table SI . Results from SPH simulations, /^j" is the radius of the donor from the stellar code in 
Rq, R[ is the largest spatial size of the pre-CEE donor in the SPH simulation after relaxation in 
Rq, a is the initial orbital separation in Rq, Pqa is the initial orbital period (in days) of the relaxed 
binary in the SPH code, /sync gives the synchronization of the red giant with the orbit (0 is for no 
rotation, 1 is for fully synchronized), A'part is the number of particles representing the giant star, 
and /rlof is the fraction of the radius of the primary with respect to its Roche lobe radius. MS 
is a main sequence star represented with 20 000 particles in mn351, and with 2000 particles 
in ms372 and ms376. PM is a point mass. The remaining columns give results and derived 
quantities: munb is the unbound mass in Mq, is the kinetic energy at infinity in 10"^^ ergs; Lp is 
the plateau luminosity in L© calculated from equation (1 ) of the main text; and tp is the duration 
of plateau phase in days, calculated from equation (2) of the main text. 
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Movie SI. Column densities in simulation ps334, as viewed perpendicular to the orbital 
plane and at times after the merger. The merger causes an outflow of material, with the field of 
ejecta dropping to lower and lower densities and column densities as it expands. At late times 
in the simulation, while the ejecta continues streaming outward, the bound merger product in 
the center of the system gradually reaches larger column densities as more material falls back 
to its surface. Visualizations generated using SPLASH (74). 

Movie S2. Temperature cross section in the orbital plane. Like visualization S1 , the scenario 
depicted is for simulation ps334, again viewed perpendicular to the orbital plane at times after 
the merger. Only temperatures between 4000 K and 5000 K are shown, corresponding roughly 
to the recombination temperature and hence roughly to the location of the photosphere. Note 
that between ~ 1 65 and 1 85 days, lines of sight from an observer looking along the orbital plane 
do not fully penetrate down to the innermost closed 4500K surface: this enlarged photosphere 
is the main cause of the a plateau phase. Near ~ 1 85 days, the shroud dissipates and the inner 
bound merger product would be revealed, corresponding to the decline stage after the plateau. 
Visualizations generated using SPLASH (74). 
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